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Sixteen strains of symbiotic bacteria from root nodules of Glycine max grown in Ottawa, Canada, 
were previously characterized and placed in a novel group within the genus Bradyrhizobium. To 
verify their taxonomic status, these strains were further characterized using a polyphasic 
approach. All strains possessed identical 1 6S rRNA gene sequences that were 99.79 % similar to 
the closest relative, Bradyrhizobium liaoningense LMG 1 8230 T . Phylogenetic analysis of 
concatenated atpD, glnll, recA, gyrB, rpoB and dnaK genes divided the 1 6 strains into three 
multilocus sequence types that were placed in a highly supported lineage distinct from named 
species of the genus Bradyrhizobium consistent with results of DNA-DNA hybridization. Based 
on analysis of symbiosis gene sequences (nodC and nifhf), all novel strains were placed in a 
phylogenetic group with five species of the genus Bradyrhizobium that nodulate soybeans. The 
combination of phenotypic characteristics from several tests including carbon and nitrogen 
source utilization and antibiotic resistance could be used to differentiate representative strains 
from recognized species of the genus Bradyrhizobium. Novel strain 0099 T elicits effective 
nodules on Glycine max, Glycine soja and Macroptilium atropurpureum, partially effective nodules 
on Desmodium canadense and Vigna unguiculata, and ineffective nodules on Amphicarpaea 
bracteata and Phaseolus vulgaris. Based on the data presented, we conclude that our strains 
represent a novel species for which the name Bradyrhizobium ottawaense sp. nov. is proposed, 
with 0099 T (=LMG 26739 T =HAMBI 3284 T ) as the type strain. The DNA G + C content is 
62.6 mol%. 



Soybean [Glycine max (L.) Merrill] is an economically 
important grain legume that can form a nitrogen-fixing 
association with species of soil bacteria belonging to the 
genus Bradyrhizobium. In a previous study (Tang et al, 
2012), populations of symbiotic bacteria associated with 
soybeans at two field sites in eastern Canada were charac- 
terized on the basis of multilocus sequence analysis (MLSA) 
of six protein encoding (housekeeping) genes. Phylogenetic 
analyses resulted in the identification of several novel 
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Abbreviations: ML, maximum-likelihood; MLSA, multilocus sequence 
analysis; ST, sequence type, RAPD, random amplified polymorphic DNA. 

The GenBank/EMBL/DDBJ accession number for the 1 6S rRNA, 
atpD, glnll, recA, gyrB, rpoB, dnaK, nodC and nifH gene sequences of 
strain 0099 T are JN 186270, HQ455212, HQ587750, HQ587287, 
HQ873179, HQ587518, JF308816, HQ587980 and JN 186287, 
respectively. Accession numbers for the 1 6 strains of Bradyrhizobium 
ottawense are given in Table S1. Strain 0085 (=LMG 26738=HAMBI 
3283) also deposited in Culture Collections. 

Six supplementary figures and six supplementary tables are available 
with the online version of this paper. 



lineages within the genus Bradyrhizobium. One of these 
novel lineages was encountered at only one of the two field 
sites and consisted of a group of 16 strains. In this work, 
we further characterize these strains using a variety of 
genotypic and phenotypic methods, and, based on the 
results, a novel species, Bradyrhizobium ottawaense sp. nov. 
is proposed. 

The 16 novel strains of symbiotic bacteria (Table SI, 
available in the online Supplementary Material) were 
obtained from root nodules of soybeans grown at a field 
site in Ottawa, Ontario (Tang et al, 2012). Bacteria were 
grown on yeast-extract mannitol (YEM) agar medium 
(Vincent, 1970) and pure cultures maintained in 20 % (w/v) 
glycerol at — 80 °C. Bacterial cells were Gram-stain-negative 
(Powers, 1995), and, based on the Schaeffer-Fulton staining 
method (Hussey & Zayaitz, 2012), were non-spore-forming. 
Colonies on YEM agar were mucoid, beige, translucent, 
circular and measured <1 mm in diameter after 7 days at 
28 °C. Cell morphology was investigated using a transmis- 
sion electron microscope (H-7000; Hitachi). Bacteria were 
cultured in stationary YEM broth for 2 days at 28 °C and 
stained with 1 % phosphotungstic acid (pH 7.0) (Hayat & 
Miller, 1990). All tested strains (0099 T , OM9 and 0085) 



065540 © 2014 Her Majesty the Queen in right of Canada as represented by the Minister of AAFC Printed in Great Britain 
3202 This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.Org/licenses/by/3.0/). 



Bradyrhizobium ottawaense sp. nov. 



had rod-shaped cells, subpolar flagella and a cell size 
(Fig. SI) that is consistent with the characteristics of the 
genus Bradyrhizobium (Garrity et al, 2005). Cell motility 
was demonstrated using the semisolid medium puncture 
method (Shields & Cathcart, 2012) with bacteria grown on 
semisolid YEM agar medium. Production of an alkaline 
reaction on YEM agar containing bromothymol blue 
after 21 days at 28 °C (Bromfield et al, 2010) and mean 
generation times (12-13 h) of bacterial cultures in YEM 
broth (Itakura et al, 2008; Wittwer, 2014) (Table 1) were 



also typical of the genus Bradyrhizobium (Garrity et al, 
2005; Vincent, 1970). 

Almost full-length 16S rRNA gene sequences were generated 
using primers 16Sa and 16Sb (van Berkum & Fuhrmann, 
2000) and sequence alignment was carried out using the 
Infernal secondary-structure-based aligner implemented in 
the Ribosomal Database Project program version 11.1 (Cole 
et al, 2014). Partial sequences of housekeeping (atpD, glnll, 
recA, gyrB, rpoB and dnaK) and symbiotic (nodC and nifH) 



Table 1. Phenotypic characteristics of Bradyrhizobium ottawaense sp. nov. strains 0099 T , OM9 and 0085 and type strains of 
species of the genus Bradyrhizobium 



Strains: 1, Bradyrhizobium ottawaense sp. nov. 0099 T ; 2, B. ottawaense sp. nov. OM9; 3, B. ottawaense sp. nov. 0085; 4, B. cytisi LMG 25866 1 ; 
5, B. yuanmingense LMG 21827 T ; 6, B. iriomotense LMG 24129 T ; 7, B. liaoningense LMG 18230 T ; 8, B. betae LMG 21987 T ; 9, B. japonieum USDA 6 T ; 
10, B. elkanii USDA 76 T . 
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*Values are based on duplicates. 
tValues are based on five replicates. 

tUninoculated control, pH 6.76 + 0.06; Ensifer meliloti ATCC 9930 T , pH 5.75 + 0.05. 
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genes were generated for the 16 novel strains of the genus 
Bradyrhizobium as well as for reference taxa not available 
in public databases. Preparation of genomic DNA, amp- 
lification, primers, nucleotide sequencing, sequence align- 
ment and editing was as described previously (Tang et al, 
2012). The dnaK sequence of 'Bradyrhizobium retamae' 
Rol9 required designing a second forward sequencing 
primer (FdnaK_SCl: GAGCAGCAGATCCGGATTCA) in 
order to obtain the 3' end of the sequence. GenBank 
accession numbers of nucleotide sequences are given in 
Tables SI and S2. 

Bayesian phylogenetic analyses were carried out using 
MrBayes version 3.2.1 with default priors (Altekar et al. 
2004). For each dataset, two concurrent analyses with 
four chains (three heated and one cold) were run for ten 
million generations with sampling every 2000 generations. 
Convergence was judged satisfactory when the average SD 
of split frequencies fell below 0.01 and the potential scale 
reduction factor statistics were approaching 1.0 (Ronquist 
et al, 2012). For each dataset, trees from the first 25% 
of sampled generations were removed as burn in and a 
majority rule consensus tree was estimated based on pooled 
post-burn-in trees. 

Best-fit substitution models were selected using the Bayesian 
information criterion implemented in jModelTest version 2 
(Darriba et al, 2012). Maximum-likelihood (ML) phylo- 
genetic analyses (Guindon et al, 2010) were carried out as 
previously described using 1000 non-parametric bootstrap 
replications to assess support (Tang et al, 2012). In all 
instances, tree topologies from Bayesian and ML analyses 
were similar; for brevity only the Bayesian trees are shown. 

Consistent with a previous report (Wang et al, 2013), two 
major groups of species of the genus Bradyrhizobium were 
evident in phylogenetic trees of 16S rRNA gene sequences: 
one represented by Bradyrhizobium japonicum and the 
other by Bradyrhizobium elkanii (Fig. S2). All 16 novel 
strains had identical 16S rRNA gene sequences and were 
placed in the phylogenetic group represented by B. 
japonicum. The sequence similarities of novel strains varied 
between 98.44 (Bradyrhizobium denitrificans) and 99.79% 
{Bradyrhizobium liaoningense) relative to type strains of the 
15 species in the B. japonicum group (Table S3). 

MLSA of at least five housekeeping genes is used as a 
reliable method to define phylogenetic relationships and to 
identify novel lineages within the genus Bradyrhizobium 
(Rivas et al, 2009; Tang et al, 2012). The Bayesian 
phylogenetic tree of six concatenated housekeeping gene 
sequences (length 3210 bp) for the 16 novel strains and 
reference taxa is shown in Fig. 1. Consistent with the results 
of Tang et al, (2012), the novel strains comprised three 
multilocus sequence types (STs) (represented by strains 
0099 T , 0085 and OM9) that were placed with high 
confidence (>95% posterior probability) in a lineage 
that was distinct from named species of the genus 
Bradyrhizobium. Relative to 0099 T , strains 0085 and 
OM9 had sequence similarities for the six concatenated 



housekeeping genes of >99.60% whereas those for 19 
named species varied between 88.35 and 95.42 % (Tables S3). 

To further analyse genetic differences between novel strains, 
we generated random amplified polymorphic DNA (RAPD) 
fingerprints for representative strains 0099 T , OM9 and 
0085 and reference taxa using four random primers (PI, P2, 
P3 and P5) and amplification methods described by Sikora 
et al. (2002). An example of the fingerprint profiles generated 
by one of the primers (PI) is shown in Fig. S3. A dendrogram 
based on the combined character matrix of fingerprint 
profiles generated by the four primers was reconstructed 
using upgma and the Dice coefficient implemented in 
GelCompare II software version 5.10 (Applied Maths). The 
three representative strains were readily distinguished and 
were placed in a single cluster separate from reference taxa 
(Fig. S4), consistent with their assignment to distinct STs, 
based on MLSA of protein encoding genes. 

DNA-DNA hybridization experiments were performed 
as described by Willems et al, (2001) using three repre- 
sentative strains (0099 T , OM9 and 0085) and relatives, 
Bradyrhizobium yuanmingense LMG 21827 T and B. 
liaoningense LMG 18230 T . DNA-DNA hybridization values 
(Table S4) among the three novel strains varied between 
73 % and 88 % whereas values for these strains hybridized 
with B. liaoningense LMG 18230 T and B. yuanmingense 
LMG 21827 T were between 23 % and 38 %. These data are 
consistent with the value of 70 % DNA-DNA relatedness 
that is considered the threshold for novel species definition 
Wayne et al, (1987). The G + C mol% content of DNA 
for the three strains, determined by HPLC (Mesbah et al, 
1989), varied between 62.3 and 63.5% which is within the 
range of DNA G + C content for the genus Bradyrhizobium. 

Analysis of accessory genes encoding symbiotic functions 
(nod and nif) can provide useful information on the host 
range and specificity of symbiotic bacteria (Laguerre et al, 
2001). The phylogenetic tree of partial nodC sequences 
of novel strains and reference taxa is shown in Fig. S5. All 
16 novel strains had identical nodC sequences and were 
placed in a group with five named species of the genus 
Bradyrhizobium that are symbiotic with soybeans. Sequence 
similarities of novel strains relative to these five species 
were 5=99.86% (Table S3). Analyses of nifH gene 
sequences (Fig. S6 and Table S3) provided almost identical 
results. 

Plant tests were carried out using Leonard jars as described 
by Tang et al, (2012). The results (Table S5) show that 
representative strain, 0099 T , elicited effective nitrogen 
fixing nodules (fix + ) on Glycine max, Glycine soja and 
Macroptilium atropurpureum, partially effective nodules 
(fix + ) on Desmodium canadense and Vigna unguiculata, 
and ineffective nodules (fix-) on Amphicarpaea bracteata 
and Phaseolus vulgaris. 

For analysis of fatty acids, strains were grown on YEM agar 
at 28 °C. Bacteria were harvested and fatty acids extracted 
as described by Sasser (1990). Fatty acid identification was 



3204 



International Journal of Systematic and Evolutionary Microbiology 64 



Bradyrhizobium ottawaense sp. nov. 



6. yuanmingense LMG 21 827 T 

6. huanghuaihaiense CCBAU 23303 T 

S. daqingense CCBAU 15774 T 
B. pachyrhizi LMG 24246 T 



B. elkanii USDA 76 T 
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Fig. 1. Bayesian phylogenetic tree of atpD-glnll-recA-gyrB-rpoB-dnaK concatenated gene sequences (3210 bp) for 
Bradyrhizobium ottawaense sp. nov. and reference taxa of the genus Bradyrhizobium (GTR + G + I substitution model). Only 
posterior probabilities ^95% are shown. Values in parentheses represent the number of strains in each of the three STs of 6. 
offawaense. Bar, expected substitutions per site. 



done using the Sherlock Microbial Identification System 
(MIDI) version 6.0 and the RTSBA6 database. Fatty acid 
profiles of three novel strains and reference taxa are shown in 
Table S6. The fatty acid profile of 0099 T was characteristic 
of the genus Bradyrhizobium with fatty acids 16 : 0 (12.83 %) 
and summed feature 8 (76.99 %) predominating (Tighe etal, 
2000). Five fatty acids (16:0, 16: 1oj5c, 18:0 and summed 
feature 8) were common to all strains tested whereas 
seven (12:0, 14:0, 17:1oj8c, 18 : lco5c, 19:0 cyclo co8c and 
summed feature 3) were detected in only some strains. 

Phenotypic characteristics were further investigated using a 
variety of tests. Utilization of carbon and nitrogen sources 
was tested using YNB and YCB (Becton Dickinson) basal 
liquid media as described by Chahboune et al, (2011) 
except that bromothymol blue was omitted. Cultures were 
incubated at 28 °C for 20 days on a rotary shaker when 
bacterial growth (turbidity) was recorded visually and 
using a Spectronic 21 UV spectrophotometer (Milton Roy) 
at 660 nm. 

Tests of intrinsic antibiotic resistance were done according 
to Bromfield et al. (2010) using YEM agar medium amended 
with filter-sterilized solutions of the following antibiotics 



(Sigma Aldrich) at final concentration (|ig ml ): erythro- 
mycin (0, 50, 100, 200), penicillin G sodium salt (0, 10, 25), 
kanamycin sulphate (0, 25, 50, 100), tetracycline (0, 5, 10, 
20), chloramphenicol (0, 25, 50), cefuroxime sodium salt 
(0, 5, 15, 30). Bacterial growth was recorded after 7 days at 
28 °C. Tests for bacterial growth on YEM agar at 10, 28 
and 37 °C, at pH 5 and pH 10 and in the presence of 1 % 
(w/v) NaCl were done as described by Xu et al, (1995). The 
combination of phenotypic characteristics listed in Table 1 
could be used to differentiate novel strains from recognized 
species of the genus Bradyrhizobium. 

Based on data for genotypic and phenotypic analyses 
presented in this study, we propose that the novel group of 
16 strains represent a novel species, named Bradyrhizobium 
ottawaense sp. nov. 

Description of Bradyrhizobium ottawaense 
sp. nov. 

Bradyrhizobium ottawaense (ot.ta.wa.en'se. N.L. neut. adj. 
ottawaense of or belonging to Ottawa, Canada). 

Cells are motile with subpolar flagella, Gram-stain-negative, 
aerobic, non-spore-forming rods (approx. 1.9 u.m long and 
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0.8 um wide). Colonies are circular, convex, translucent, 
beige, < 1 mm diameter after 7 days at 28 °C on YEM agar 
medium. Mean generation time approx. 12-13 h in YEM 
broth at 28 °C. Produce an alkaline reaction on YEM agar. 
Growth occurs at pH 5-10 (optimum, pH 7.0). Growth is 
optimal at 28 °C but no growth occurs at 10 or 37 °C or 
in the presence of 1 % (w/v) NaCl. Utilizes D-mannose, 
D-xylose, L-arabinose, D-glucose, D-mannitol, maltose, L- 
rhamnose, L-proline, L-glutamic acid and L-threonine but 
not lactose, trehalose, succinic acid, L-sorbose, aspartic acid, 
glycine or DL-a-alanine. Resistance (ug ml -1 ) to antibiotics 
for representative strain 0099 T : kanamycin (100), tetra- 
cycline (20), chloramphenicol (100), erythromycin (100), 
cefuroxime (<5), penicillin (<10). Summed feature 8 
(18 : li»6cand/or 18 Aw7c) and 16 : 0 are predominant fatty 
acids. Strain 0099 T elicits effective nodules on Glycine max, 
Glycine soja and Macroptilium atropurpureum, partially 
effective nodules on Desmodium canadense and Vigna 
unguiculata, and ineffective nodules on Amphicarpaea 
bracteata and Phaseolus vulgaris. 

The type strain, 0099 T (=LMG 26739 T =HAMBI 3284 T ), 
was isolated from an effective nodule of Glycine max in 
Ottawa, Ontario, Canada. The DNA G + C content of the 
type strain is 62.6 mol%. 
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